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Abstract. We study the formation and growth of equilibrium clusters in 
a suspension of weakly-charged colloidal particles and small non-adsorbing 
polymers. The effective potential is characterised by a short-range attraction 
and a long-range repulsion. The size, shape and local structure of the clusters are 
studied using three-dimensional particle microscopy. We observe a rapid growth in 
the mean cluster size and the average number of nearest neighbours approaching 
the gel boundary. 
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1. Introduction 

The van der Waals theory is the foundation for our current understanding of fluid 
structure and phase equilibria pQ. It predicts, for an interaction potential consisting 
of a short-range repulsion and a long-range attraction, a first-order transition (below 
T c ) between a dilute vapour and a dense fluid or liquid. Here we consider the unusual 
situation which pertains when the range of the repulsive and attractive forces are 
interchanged. 

Experiments H E] , theory El HHj and simulation [111112111 suggest 
that a combination of a very short-ranged attraction and a much longer range repulsion 
leads to a highly nontrivial phase behaviour. At low densities, there is compelling 
evidence for the existence of a fluid phase of finite-sized equilibrium aggregates or 
clusters [2 E3 ■ The size and shape of these clusters are expected ^21 to be very 
sensitive to the relative range of the attractive and repulsive components (Vn(r)) of 
the total interparticle potential. Estimates of the ground state energies of clusters of 
different sizes suggest that the equilibrium shape changes from spherical to linear as 
the cluster diameter exceeds the decay length of Vn(r) In addition recent 

experiments |H] have revealed that, at high densities, the competition inherent in 
this system can also stabilize new particle architectures. Three-dimensional confocal 
microscopy in a system of weakly charged colloidal spheres show, at high densities, 
the formation of a macroscopically percolating gel composed of linear aggregates of 
face-sharing tetrahedra of particles. The colloidal particles are arranged in the motif 
of a helical spiral in which each particle is connected to six neighbours, analogous to 
the structure first proposed by Bernal in his classic study of hard-sphere glasses |14) . 

In this paper, we present an experimental study of a carefully-characterized model 
system of colloidal spheres. The particles interact through an effective potential which 
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Figure 1. Behaviour of a charged colloid-polymer mixture. Circles represent 
samples that form isolated clusters. Diamonds depict macroscopic percolating 
gels. The solid line identifies the experimentally-determined gel boundary (f> g . 
Data from circled samples A— G are shown in figure 



consists of a short-range attraction, generated by the exclusion of a non-adsorbing 
polymer, and a long-range screened electrostatic repulsion. We focus on what happens 
to the size and shape of the self-assembled clusters as the gel boundary is approached. 

2. Experimental details 

We use a well-studied model system: random-coil polystyrene (radius of gyration 
r g = 92 nm) mixed with PMMA spheres (diameter a — 1700 nm) dispersed in 
a near-density and near-index matched mixture of 78:22 % by weight cycloheptyl 
bromide and cis-decalin. The PMMA spheres were labelled with the dye DilCis 
to make them visible by fluorescence microscopy. Exclusion of the non-adsorbing 
polymer from a region between two neighbouring particles generates a short-range 
attraction. The range and strength of which is controlled by r g and the free 
polymer concentration (f>*. For the cases studied here the width of the attractive 
range is roughly O.lcr. Measurements of the particle electrophoretic mobility using 
phase analysis light scattering revealed that the PMMA particles develop a small 
positive charge, Q « +145e, in the density-matched solvent used in this study. The 
Debye screening length was estimated from measurements of the solvent conductivity 
(170 ± 30 pS cm" 1 ) as of the order of na « 0.5. 

Mixtures with different colloid (4> c ) and polymer concentrations (0*) were 
prepared and homogenized by extensive tumbling. To minimize decomposition and to 
ensure reproducibility all samples, once prepared, were stored in light-tight containers 
at 4°C, prior to measurement. The size and internal structure of clusters were studied 
using laser scanning confocal microscopy. The suspensions were sealed in a cylindrical 
cell of ~ 50 /iL volume. Three-dimensional confocal microscopy was used to distinguish 
between a phase of isolated clusters and the formation of a macroscopic percolating 
network of particles. Our observations are summarized in figure ^ To follow the 
transition between isolated clusters and a gel phase in detail we focus here on the 
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sample sequence A-G, figure^ at 4>* p — 0.95. A stack of typically 345 images, spaced 
by 0.16 /jm vertically, was collected from a representative volume of 73 /im x 73 
(im x 55 jum within each suspension. Each stack contained, depending on (f> c , the 
images of between 1500 and 45000 spheres. The images were analyzed quantitatively 
using algorithms similar to those described in \H>\ to extract the three-dimensional 
coordinates of each sphere. The position of each sphere was measured with a precision 
of approximately 70 nm. 



3. Data analysis 

The number, size and shape of the clusters were identified from the coordinate list. We 
define particles which are 'bonded' to each other by their separation r. If r < r , where 
the cut-off distance ro was identified with the position of the first minimum in the pair 
distribution function g(r), then the particles were considered neighbours. Standard 
algorithms were used to partition particles into clusters of size s and to evaluate the 
cluster size distribution n s . The second moment of the cluster distribution (S2) was 
used to characterize the average cluster size JB] , 

(82) = ^f^. (1) 

The geometric properties of the clusters were evaluated from the eigenvalues Ax < 
A2 < A3 of the radius of gyration tensor S, defined as 

1 N 

S af 3 = ^ a ~ r iA[ r i,H ~ r jA ( 2 ) 

with J", the position of the ith sphere in the cluster of size N and (a, 0) — 1, 2, 
3 denoting the corresponding Cartesian components. We characterize the size and 
shape of the cluster in terms of the rotational invariants of the tensor S. The trace 
Y]j Aj gives the usual isotropic squared radius of gyration r 2 while the anisotropy of 
the cluster is measured by the asphericity A2 |17| . 

, _ (A 1 -A 2 ) 2 + (A 2 -A 3 ) 2 + (A 3 -A 1 ) 2 

M ~ 2(A 1+ A 2 + A 3 ) 2 • (3) 

The asphericity is normalized so that A2 ranges from zero for a spherically symmetric 
cluster to one for a rod-like cluster. 



4. Results 



The left panel of figure [21 shows the evolution in the proportion of single unbonded 
particles / as a function of <p c . The boundary between a non-percolating cluster 
phase and the percolating gel occurs at <j) g ~ 0.1. The most striking feature is the 
abrupt change in the number of single particles approaching the gel boundary. At 
low colloid concentrations / is approximately constant at / sa 0.22, independent of 
cj) c - the remaining 78% of particles being incorporated into clusters. The proportion 
of single particles however drops rapidly in the gel phase, although it remains finite 
even for dense gels (/ « 10~ 2 at <p c = 0.3, for instance). The coexistence evident 
in these samples between free and bonded particles highlights the dynamic nature of 
the assembly process. Visual observation revealed particles continuously joining and 
leaving clusters. To provide further evidence that equilibrium had been reached we 
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Figure 2. (a) The proportion of unbonded particles / and (b) the average cluster 
size (S2), as a function of colloid volume fraction C , for samples marked A— G in 
figure Q 
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Figure 3. Size dependence of the cluster radius of gyration: sample A, <f>c = 0.011 
(triangles), sample B, <j> c = 0.076 (inverted triangles), and sample C, <f> c = 0.086 
(circles). The lines are power law fits to the data. 



studied the time-dependence of /. We found no qualitative change on the period of a 
week, during which time particle sedimentation was negligible. 

Figure |2Ib) shows the growth in the average cluster size (sa), defined in Q, 
with increasing 4> c for all non-percolating states. Clearly the the average cluster size 
increases rapidly approaching the gel boundary. This rapid growth is mirrored in the 
number of particles s max contained in the largest observed cluster, which increases 
from 37 (<p c = 0.011), to 155 (cf> c = 0.076), to reach s max = 2340 at (f> c = 0.086. 
To quantify these changes in greater detail we determined the size dependence of the 
ensemble-averaged radius of gyration (r g ). For fractal aggregates r g displays a power 
law dependence on the number of particles s, (r g ) ~ s 1 ^ d/ , where df is the fractal 
dimension. Figure |3] shows plots of r g versus s, at three values of <p c . Interestingly we 
find that at all densities the clusters show fractal scaling but that the fractal dimension 
df shows a marked dependence on cf> c . At low cf> c , the fractal dimension found of 
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df — 1.73 ± 0.07 is consistent with the diffusion-limited cluster-cluster aggregation 
(DLCA) value in three dimensions (df = 1.78 |181 119j 1. While at high <fi c , the value 
of df = 2.14 ± 0.03 is consistent with the reaction-limited cluster-cluster aggregation 
(RLCA) value in three dimensions (df = 2.1 f2H). This apparent crossover between 
DLCA and RLCA kinetics with (f> c is unexpected. However it is unclear if fractal 
models for non-equilibrium growth are applicable to the equilibrium case of assembly 
studied here for which, as we indicate below, the cluster shape is not truly self-similar. 

A simple indicator of the anisotropy of the clusters is provided by the asphericity 
(^2), defined in ©. In contrast to the r g data we find little (f> c dependence of the 
average cluster asphericity, plotted in figure EJa). Both low and high density samples 
show a similar monotonic decrease in (A2) with s, consistent with a progressive change 
in shape with increasing cluster size. The clusters formed are therefore not self- 
similar. Clusters smaller than about 30-50 monomers show a definite tendency to 
be elongated in one dimension as evidenced by the increasing value of (^2) found at 
small s. This preference for a linear or rod- like shape seems to be more pronounced the 
smaller the cluster. These observations are supported by images of individual clusters 
which reveal characteristic one-dimensional bundles of particles. By contrast, while 
large clusters remain anisotropic the levels of asymmetries are significantly lower and 
typical of clusters formed from either random percolation ((A2) = 0.25 [21]) or DLCA 
models (A 2 ) = 0.32 This crossover suggests that in large clusters the short one- 

dimensional bundles have branched many times generating clusters with geometries 
controlled by random percolation. 

The tendency seen for one-dimensional aggregation at short scales might indicate 
that the cluster is constructed locally from segments of Bernal spiral structures. To 
check for this possibility, we calculated the average number of nearest neighbours (n). 
As shown in the right hand panel of figure 0J (n) grows abruptly with increasing 
4> c approaching the gel boundary. Above (j) g , the mean coordination number is 
insensitive to <j) c and is approximately six, consistent with the characteristic six-fold 
coordination of a Bernal spiral. Moreover figure 01 reveals that (n) is less than six 
in the non-percolating cluster samples (A-C). This is particularly evident in sample 
B (4> c = 0.076) which while it contains clusters of upto 100 monomers has a mean 
coordination number of only (n) — 2.7. The lack of spiral ordering, shown by the low 
values of (n), is also supported by direct inspection which reveals little evidence for 
Bernal ordering in any of the non-percolating cluster samples. 

5. Conclusions 

In this work we have explored the change in the size and shape of the clusters formed 
in a colloidal system of short range attractions and long range repulsions, approaching 
the gelation transition. In contrast to a purely attractive system the clusters are 
equilibrium entities and their properties time- independent. Our most noteworthy 
finding is that the average size (S2) of the cluster distribution, the proportion of 
unattached particles, and the mean coordination number all change abruptly as 
<p — > 4> g . Although we have studied only a limited number of concentrations (at 
one specific choice of potential parameters) these observation suggest that, for this 
system, gelation has many of the hallmarks of an equilibrium transition. Qualitatively 
our observations appear consistent with the suggestion by Groenewold and Kegel [71IH] 
that under the right conditions a fluid of isolated clusters will undergo a spinodal 
decomposition, driven by counter-ion condensation, to form a macroscopic gel. This 
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Figure 4. (a) Cluster asphericity as a function of cluster size. The dotted lines 
depict the asphericity expected for random percolation )2 1 1 and DLCA models 
1221 . (b) Average number of neighbours as a function of volume fraction rp c . 
Symbols have the same meaning as figure 131 



raises the enticing possibility of an ideal thermally-reversible gel - one in which the 
gel-line might be approached infinitely closely and the dynamics of arrest studied 
precisely and reproducibly, free from the dynamical complications of purely attractive 
systems. Further experiments are planned to test these ideas. 
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